Dynamo mechanism: Effects of correlations and viscosities 
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, We analyze the effects of the background velocity and the initial magnetic field correlations, and 

. viscosities on the turbulent dynamo and the a-effect. We calculate the a-coefficients for arbitrary 

magnetic and fluid viscosities, background velocity and the initial magnetic field correlations. We 
explicitly demonstrate that the general features of the initial growth and late-time saturation of the 
■ magnetic fields due to the non-linear feedback are qualitatively independent of these correlations. 

We also examine the hydrodynamic limit of the magnetic field growth in a renormalization group 
framework and discuss the possibilities of suppression of the dynamo growth below a critical rotation. 
We demonstrate that for Kolmogorov- (K41) type of spectra the Ekman number M ~ 1/2 for 
dynamo growth to occur. 
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; I. INTRODUCTION 

Magnetic fields are ubiquitous. All astrophysical objects are known to have magnetic fields of diS^erent magni- 
^ tudes,e.g., 1 gauss at the stellar scale to 10~^ gauss at the galactic scale [1]. The origin of such fields {primordial 



(N 



field) is not very clear - there are several competing theories which attempt to describe this [2]. However, a finite 
' magnetic field in any physical system undergoes a temporal decay due to the finite conductivity of the medium. 
(— I So, for steady magnetic fields to occur in astrophysical bodies, there has to be a mechanism of regeneration of the 
O magnetic fields, which takes place due to the dynamo process [1,3]. Most astrophysical bodies are thought to have fast 
, dynamo operating within themselves (there are exception to this, e.g., the Moon, Venus and Mars in our solar system) 

resulting into exponential growth of the magnetic fields. This mechanism requires a turbulent velocity background [1] 
I , [though non-turbulent velocity fields too can make a seed (initial) magnetic field to grow (for details see [3]), we will 
' not consider such cases here]. Since the dynamo equation, in the linear approximation (see below) gives unbounded 
exponentially growing solutions for the long wavelength (large scale) part of the magnetic fields, it is linearly unstable 
1^ in the low wavenumber limit. However, one does not see a perpetual growth of magnetic fields in the core of the earth 
^SJ or in the sun. For example, geomagnetic fields (~ 1 gauss) are known to be stable for about 10^ years [1]. Thus, the 
' physically realisable solutions of the dynamo equations cannot be unstable in the long time limit. It is now believed 
xl , that the non-linear feedback due the Lorentz force term in the Navier-Stokes equation is responsible for the saturation 
of the magnetic field growth (see, e.g., [1]). 

The study of this problem has already been the subject of previous work by many groups. For example Pouquet, 
Frisch and Leorat [4] studied the connections between the dynamo process and the inverse cascade of magnetic and 
kinetic energies within a eddy damped quasi-normal Markovian approximation. Moffatt [5] has examined the back 
'■^ ' reactions due to the Lorentz force for magnetic Prandtl number Pm, ^ 1 by linearising the equations of motion 
, of three-dimensional (3c?) magnetohydrodynamics (MHD). Vainshtein and Cattaneo [6] discussed several nonlinear 
restrictions on the generations of magnetic fields. Field et al [7] discussed nonlinear a-effects within a two-scale 
approach. Rogachevskii and Kleeorin [8] studied the effects of an anisotropic background turbulence on the dynamo 
process. Brandenburg examined non-linear a-effects in numerical simulation of helical MHD turbulence [9]. In 
particular, he examined the dependences of dynamo growth and the saturation field on the magnetic Prandtl number 
' Pm (the ratio of the magnetic- to the kinetic- viscosities). Bhattacharjee and Yuan [10] studied the problem in a 
9^ 1 two-scale approach by linearising the equations of motion. 

Dynamo mechanism has two competing processes at work: amplification of the magnetic field by the dynamo 
process and ohmic dissipation due to finite resistivity of the medium concerned. Which one among these two effects 
will dominate depends on the case in study. In some specific models, however, one can analyze this completely. A 
good example of such models is the Kraichnan-Kazantzev dynamo [11,12] where the velicity field is assumed to be 
Gaussian-distributed, delta-correlated in time and the magnetic field is governed by the Induction equation [22]. In 
this model the statistics of the velocity field is taken to be parity invariant so that the a-effect is ruled out. The 
main results from this model include i)the existence of dynamo in the infinite magnetic Reynolds number limit for a 



o 

e 

I 



o 
o 



X 



1 



particular choice of the variance of the velocity distribution [13] and ii)the existence of a critical magnetic Reynolds 
number only above which dynamo growth is possible [14]. However, not much is known about this when invariance 
due to parity is broken and when the velocity field is not temporally delta-correlated. In a recent simulations [15] the 
authors found, in a model simulation for the solar convection zone, a monotonic increase of the horizontal a-effect 
with rotation. Kida et al showed, in numerical simulations, that unless magnetic hyperviscosity is less than a critical 
value, magnetic fields did not grow [31], confirming the existence of a critical magnetic Reynolds number (Rm)- 

Our studies generalize the existing results. In this paper we use a minimal model of a-cfFect (see below) to study 
dynamo with a-effect to calculate the a coefficient for arbitrary correlations and viscosities, and ask the following 
questions: 

1. Do the turbulent dynamo growth and the saturation processes require any turbulent background? Or do they 
function with arbitrary parity-breaking and fluctuating velocity and initial magnetic field correlations? ^ . 

2. What is the hydrodynamic limit (long wavelength limit) of the dynamo problem? By this we ask how the 
magnetic field correlations scale in the infra red limit during the initial-growth regime. 

3. Can arbitrarily large magnetic viscosity prevent dynamo growth? In other words, is there a critical magnetic 
Reynolds number Rm above which the dynamo growth sets in? 

To study the above mentioned questions we employ a diagrammatic perturbation theory, which has been highly 
successful in the contexts of critical dynamics [18], driven systems [19], etc. This can be easily extended to higher 
orders in perturbation expansion and is very suitable for handling continuous kinetic and magnetic spectra, unlike the 
two-scale approximation. This was first used to study stationary, homogeneous and isotropic MHD in Ref. [20]. We use 
this method to study non-stationary statistical states (dynamo growth) which facilitates studies on the hydrodynamic 
limit of the dynamo problem in a renormalisation group framework. We use diagrammatic perturbation theory to 
calculate expressions for the a coefficiants for arbitrary background velocity and initial magnetic field correlations 
and magnetic Prandtl number Pm for both the early growth and the late time saturation. With our expressions for 
a we examine the three issues mentioned above. 

We investigate these for arbitrary correlations and magnetic Prandtl number with no approximations other 
than the existence a perturbation theory. Our principal results are: 

• We calculate the a-coefficients for arbitrary correlations and viscosities. 

• We examine the hydrodynamic limit in the kinematic regime and predict the existence of a critical Rm or 
rotation above which dynamo growth will occur for certain correlations with infra red singularity. 

In our all our studies, we do not assume any variance for the velocity field. Instead, we use the Navier-Stokes equation 
to describe the dynamics of the velocity field. This allows us to use a renormalisation group framework to study the 
hydrodynamic limit. 

The first question that we investigate is phenomenologically very important because different systems may have 
different velocity and initial magnetic field spectra. Therefore, it is important to imdcrstand the dependence of the 
dynamo on these spectra. We explicitly demonstrate that the nonlinear feedback of the magnetic fields on the velocity 
fields in the form of the Lorentz force stabilises the growth for arbitrary velocity and initial magnetic field correlations. 
This demonstrates that the basic features of the dynamo mechanism arc qualitatively independent of the velocity and 
magnetic field spectra and, essentially, are properties of the SdMHD equations. Details (e.g., the values of the a- 
coefficients) of course, depend upon the actual forms of the spectra. Our renormalization group analysis indicates that 
dynamo growth takes place only if the Ekman number Af < 1/2 (for a given Rm) when the velocity and the initial 
magnetic field spectra are sufficiently singular in the long wavelength limit. The structure of this paper is as follows: 
In Sec. II we discuss the general dynamo mechanism within the standard linear approximation for arbitrary velocity 
and initial magnetic field correlations and viscosities. In Sec. Ill B we show that beyond the linear approximation 
non-linear effects lead to the eventual saturation of magnetic field growth for arbitrary background kinetic energy 
and initial magnetic energy spectra, and viscosities. We elucidate how different background kinetic energy and initial 



^By a turbulent background we do not mean any kind of fluctuating state but a fluctuating state with Kolmogorov (K41) 
spectra oc k~^^^ for the kinetic and magnetic energies and cascades of appropriate quantities; if there is no mean magnetic fleld 
then the energy spectra is expected to be K41-type - see Ref. [16]. 
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magnetic energy spectra affect the values of the a-coefficients. In Sec. IV we analyze the initial dynamo growth in 
a renormalization group framework. We show that for sufficiently singular velocity and magnetic field spectra the 
Ekman number must be < 1/2 for the magnetic fields to grow. For velocity and magnetic field spectra which go to 
zero in the long wavelength limit there are no such restrictions. In Sec.V we present our conclusions. 



II. DYNAMO GROWTH: THE LINEAR APPROXIMATION 



In the kinematic approximation [1,21], i.e., in the early-time regime, when the magnetic energy is much smaller than 
the kinetic energy (J u^d^r >> J b^d^r, where u(r, t) and b(r, t) are the velocity and magnetic fields respectively) the 
Lorentz force term of the Navier Stokes equation is neglected. In that weak magnetic field limit, which is reasonable 
at an early time, the time evolution problem for the magnetic fields is a linear problem as the Induction equation [22] 
is linear in magnetic fields b: 

^ = Vx(uxb)+/xV2b, (1) 

ot 

where /i is the magnetic viscosity. The velocity field is governed by the Navier-Stokes equation [23] (in the absence of 
the Lorentz force) 

^ + u.Vu = -^+i/V2u + f. (2) 
at p 

Here v is the fluid viscosity, f an external forcing function, p the pressure and p the density of the fluid. We take f 
to be a zero mean, Gaussian stochastic force with a specified variance (see below). 

In a two-scale [1] approach one can then write an effective equation for B, the long- wavelength part of the magnetic 
fields [1]: 

f)TK 

— = Vx(UxB) + VxE + /xV2B, (3) 

ot 

where the Electromotive force E = (u x b). U is the large scale component of the velocity field u. An Operator 
Product Expansion (OPE) is shown to hold [21] which provides a gradient expansion in terms of B for the product 
E = (u X b) [1] 

dB 

Ei=a^jBj+(3ijk-g^ + .... (4) 

For homogenous and isotropic flows {aij = aSij) Eq.(4) gives, 

/)B 

— =Vx(UxB) + aVxB + /xV^B, (5) 

ot 

which is the standard turbulent dynamo equation. Here p, now is the effective magnetic viscosity which includes both 
the microscopic magnetic viscosity and the turbulent diffusion, represented by Pijk in Eq.(4). a depends upon the 
statistics of the velocity field (or, equivalently, the correlations of f). Retaining only the a -term and dropping all 
others from the RHS of Eq.(5), the equations for the cartesian components of B become (we neglect the dissipative 
terms proportional to as we are interested only in the long wavelength properties) 




B^{Kt)\ ( -k, ky \ /B,{k,t)\ 

By(k,t) \ =ia\ k, -fc^ By{k,t) 
B,(k,t)J \-ky k^ J \B,ik,t)J 



The eigenvalues of the matrix is A = ±ik, 0. Thus depending on the sign of the product ak, one mode grows and the 
other decays. The third mode is unphysical, because the corresponding eigenfunction is proportional to k and hence 
in confiict with V • B = 0. Since growth rate is proportional to and dissipation is proportional to fc^, large scale 
fields continue to grow leading to long wavelength instability. Thus in the long time limit effectively only the growing 
mode remains. Growth rate a is a pseudo-scalar quantity, i.e., under parity transformation r — r, a —a [1,21]. 
Since a depends upon the statistical properties of the velocity field, its statistics should not be parity invariant. This 
can happen in a rotating frame, where the angular velocity explicitly breaks reflection invariance. 
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III. FORMULATION OF THE DYNAMO PROBLEM IN A ROTATING FRAME 



The Navier-Stokes (NS) (including the Lorentz force) and the Induction equation in an inertial frame in (k, t) space 
take the form 

+ ip,,p(k) J2 t)M^ - q, = IPmi^) ^)^3'(^ - q- i) - ^k^u^ + 0, (6) 

9 Q 

and 

^h!^ = .^(k) ^ „^.(q, t)bp{k - q, i) - nk%i. (7) 

Q 

Here, Ui{k,t) and bi{k,t) are the fourier transforms of Ui{r,t) and bi{r,t) respectively, Pijp(k) = Pij{k)kp + 
Pip{k)kj, Pijp{k) = Pij{k)kp — Pip{k)kj, Pij is the projection operator, which appears due to the divergence-free 
conditions on the velocity and magnetic fields (we consider incompressible fluid for simplicity). Equations (6) and 
(7) have to be supplemented by appropriate correlations of fi and initial conditions on 6j. We choose fi{k,t) and 
bi{k, t = 0) to have zero mean and to be Gaussian distributed with the following variances: 

(/,(k, t)fj{-k, 0)) = 2PMk)6{t), (8) 

{h{k,t = 0)bj{-k,t = 0)) = 2PijD2{k), (9) 

where Di and D2 are some functions of k (to be specified later). 

In a rotating frame with a rotation velocity = Oz the Eq.(6) takes the form 

+ 2{n X u)i + ^Pijp{k) ^ Uj{q, t)up{k -q,t) = ^P^Jp{k) ^ 6, (q, t)bp{k - q,t) + vV'^u, + /^(k, t), (10) 

9 9 

whereas Eq.(7) has the same form in the rotating frame, x u is the coriolis force. The centrifugal force x (fi x r) 

is a part of the effective pressure=p + x which does not contribute to the dynamics of incompressible flows. 
The bare propagator G„ (obtained from the linearized version of Eq.(lO)) of Ui 

iu)+vk^ 20 n 

Q — I 20 iu+uk^ rv 

1 

such that u = Gu f where u is the column vector 




u 



One can verify that with the form of the bare propagator given above, an odd-parity part in the velocity auto-correlator 
{ui{k,t)uj{—k,Q)) appears which is proprotional to the rotation O. Notice that G"' is different from GJ^'^^ - this 
is just the consequence of the fact that distinguishes the z-direction as a preferred direction in space, making the 
system anisotropic. However for frequencies lo » 57 or length scales fc^ >> O (here z is the dynamical exponent) 
isotropy is restored. In that regime, to 0{Vt) the role of the global rotation is to introduce a non-zero odd-parity part 
in {uiUj) proportional to 17. This can be also seen by noting that in the inertial frame the correlation {ui(k)uj{—k)) 
is of the form Pij{k)A{k) [cf. Eq.(8)] where A{k) is a scalar function of k and hence in the rotating frame the 
correlator is proportional to RPijBF where Ra are appropriate rotation matrices (we have suppressed the indices). 
Similarly, initial magnetic field correlations, given by Eq.(9) transforms accordingly in the rotating frame. Since 
rotation matrices act on {ui(k)uj{—k)) and Eq.(9) in the same way, magnetic field auto-correlator {bi{k,t)bj{—k,0)) 
has an odd parity part in the rotating frame with the same sign as the odd parity part in the velocity correlator. 
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Thus the effects of rotation can be modeled (to the lowest order) by introducing parity breaking parts in Eqs.(8) and 

(9) [1] 

(/,(k,t)/,(-k,0)) = 2P,jDi{k)Sit) + 2ie,,pkpDi{k)6{t), 
{bi(k, t = 0)6j(-k, i = 0)) = 2PijD2{k)S{t) + 2i€ijpkpD2{k), (11) 

in conjunction with the Eqs.(6) and (7), where e^p is the totally antisymmetric tensor in 3d. This way of modeling 
rotation effects is, of course, only approximate, but suffices for our purposes as this explicitly incorporates parity 
breaking. One can, however, construct experimental set ups [1] which are described correctly by Eqs.(ll). The 
parity breaking parts in the noise correlations or initial conditions ensure that the velocity and the initial magnetic 
field correlators have non-zero odd parity parts, as would happen in a rotating frame. An important dimensionless 
number is the Ekman number M = which can be related to D\ by equating the parity braking parts of the velocity 
correlator calculated from (linearized) Eq.(lO) and Eq.(8) with that from Eqs. (6) and (11). This gives Di = 2M~^Di. 
Now, one may ask what is the relative sign between Di and 1)2? Since the parity breaking parts of the correlators of 
the velocity and the magnetic fields have same sign and arc proportional to D\ and £>2 respectively, D\ and D2 must 
have same sign. As already noted, introduction of parity breaking terms in the force/initial correlations is well-known 
in the literature, we, nevertheless, give the analysis in details in order to emphasise on the fact that fiuid and magnetic 
helicities must have the same sign. Furthermore, for a complete description of the effects of rotation, in addition to 
the coriolis force, a forcing with a preferred direction is also required. We, however, do not include all these details 
as introduction of parity-breaking correlations is sufficient for our purposes. In this sense, this can be thought of as 
a reduced or a minimal model for dynamo. One may note that a nonzero kinetic helicity is required for the a-effect 
as the a-cocfficicnt is proportional to the kinetic helocity. Even though a global rotation explicity breaks the parity 
invariance of the system under space reversal, rotation alone is not enough to yield a non-zero helicity. This is because 
the helicity is pseudo-scalar and, therefore, can be constructed only out of an axial vector (here, rotation n) and a 
polar vector. In typical astrophysical settings, the latter one could be provided by, say, a density inhomogeneity. Even 
though this is not contained in Eq. (6), our minimal model, nevertheless, produces a finite helicity due to the helical 
nature of the forcing function. Thus, our minimal model is able to capture both the breakdown of parity due to the 
rotation and the generation of helicity due to the rotation and any other preferred direction. 

A. The a in the kinematic approximation: Dependences on background velocity and initial magnetic field 

spectra 

In the kinematic approximation, which neglects the Lorentz force term of the Navier-Stokes equation, the time 
evolution of the magnetic fields follows from the linear Induction Equation (1). We assume, for the convenience of 
calculations, that the velocity field (u) has reached a statistical steady state. This is acceptable as long as the loss 
due to the transfer of kinetic energy to the magnetic modes by the dynamo process is compensated by the external 
drive. In the kinematic (i.e., linear) approximation, we work with the Eqs. (6) (without the Lorentz force) and (7). 
We choose /;(k, t) to be a zero-mean, Gaussian random field with correlations 

(/,(k,f)/„,(k,0)) = 2Pi^D^{k)5{t)+2iei^nknDi{k)6{t). (12) 

Our initial conditions for the magnetic fields are 

(6a(k, t = 0)6/3(-k, t = Q)) = 2P^pD2{k) + 2iea0^k^D2{k), (13) 

Since we arc interested to investigate the dynamo process with arbitrary statistics for the velocity and magnetic 
fields we work with arbitrary Di{k), I)i{k), D2{k) and D2{k). For K41-type spectra, we require [24] Di{k) = 
Dik~^,Di{k) = Dik~'^,D2{k) = D2k~^/^ and D2{k) = k~^^^. These choices ensure that under spatial rescaling 
X ^ /x, u, b ^ l^/^{u, b} which is the Kolmogorov scaling [24]. Note that both the force correlations in the Eq.(6) 
and the initial conditions on Eq.(7) have parts that are parity breaking, in conformity with our previous discussions. 
We now calculate the a-term. We use an iterative perturbative method which is very similar to and discussed in 
details in Ref. [19]. In this method, terms in each order of the perturbation series can be represented by appropriate 
Feynman diagrams [19]. Even though, for simplicity, we confine ourselves to the lowest order in the perturbation 
theory (represented by the tree level diagrams), which is sufficient for our purposes, higher order calculations rep- 
resented by higher order digrams can be done in a straight forward manner. Below we give the expression for a 
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in the kinematic approximation (which we caU the 'direct' term - responsible for growth) in the lowest order of the 
perturbation theory (see Fig. la): 

((u X h)^)D = { / (fi0-fUi3{ci, i)6T,(k - q, t)) = ( / eafjjU/sici, t)ejsxi{^ - q)5'W,,(qi, ii)&r (k - q - qi, ii)Go(fc -q,t- h)) 

Jq,qi Jq 

(14) 



from which one can read the a-term: 



(k,i)= / ^^^2^^ e0npqpea0'r{-i)q5bpiKt = 0)[ 

Jq 



1 _^exp(-2tV)i 



q^{v + li) q^{v-ii) 



(15) 



giving an = ^f^ v{v+ij,) Iq '^llT+^q^ ^'^^ large t. The sufSx D refers to growth or the direct term, as opposed to feedback 
which we discuss in the next Sec.IIIB. The growth term is proportional to \k\ and diffusive decay proportional to fc^. 
The angular brackets represent averaging over the noise and initial- condition ensembles. 



D, 



o 



k-q k k-q k 

FIG. 1. Tree level diagrams for < u(q) x b(k — q) >. (a)Contribution to growth term ao'- A solid line indicates a bare 
magnetic field response function, a broken line indicates a bare velocity response function, a 'o' joined by two broken lines 
indicates a bare velocity correlation function (proportional to Di), a wavy line indicates a magnetic field, a solid triangle 
indicates a ub vertex, (b) Contribution to feedback term ap: A solid line indicates a bare magnetic field response function, a 
broken line indicates a bare velocity response funtion, a 'o' joined by two broken lines indicates a bare magnetic field correlation 
function (proportional to i32), a wavy line indicates a magnetic field, a solid triangle indicates a ub vertex. 
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B. Suppression of growth rate: Nonlinear feedback 



When the magnetic fields become strong, it is no longer justified to neglect the feedback of the magnetic fields in 
the form of the Lorentz force. So we need to work with the full Eqs.(6) and (7). The ideas of OPE as elucidated 
in Sec. II are still valid for the full non-linear problem. But the value of a is expected to change from its value in 
the linear problem. In presence of the Lorentz force there is an additional contribution to a (Fig. lb). To evaluate 
that, we follow a diagrammatic perturbation approach similar to that described in the previous Section. Here also we 
restrict ourselves to the lowest order only (i.e., the tree level diagrams) though extension to higher orders is straight 
forward. We obtain 

((u X bi)^) = ( / eijpUj{q,t)bp{\s.- q,t)) (16) 

Jq 

I Pjmn{q)Gl{q,t-ti) 6„(qi,ti)6„(q-qi,ii)C?^(fc-g,t)6p(k-q,i = 0)) (17) 

Jq 



which gives {F refers to feedback) 

Jq 2aD\q\ - ^nq-" 

which, after some simplifications, yields, 

r,2 



(18) 



2^3 4 f D2{q,t)q' 

apit) = —— / T- — 2' (19) 

3 15 Jg aD\q\ - 2jj,q^ 

where D2{q,t) = exp[2aD\q\t — 2iiq^t]D2{q) is a growing function of time for small wavenumbers. As before, angular 
brackets refer to averaging over noise and initial-condition ensembles. Thus ap grows in time. 

Since, at any finite time i, when the non-linear feedback on the velocity field due to the Lorentz force is nolonger 
negligible, both an and ap are non-zero and we get 

<Pq Di{q) 



2S3 f d-'q 

- ^± f ^ 

3 15] (27r 



[\{aD+ap)q\ - {u + ij.)q'^]'- 



d^q D2{q,t)q 



2 



(20) 



0^ \{<^D + ap)q\ -2/ig2' 
with 

D2{q,t) = exp[2(az, - ap{t))\q\t - 2fiqH]D2iq). (21) 

Equations (20) and (21) are to be solved self-consistently [17]. Thus the net growth rate is proportional to \(aD + ap)k\ 
for the mode i?j(k, f). The expressions (20) have apparent divergences at finite q; so in perturbative calculations one 
should treat the a-tcrms as perturbations which remove these divergences. This problem is akin to that in Kuramoto- 
Shivashinsky equation for flame front propagation [25]. So the expressions for ao and ap are 



2_^ r dhj Dijq) 
3 J {2nr i^iiy + fi)q^ 



253 4 f d^q D2{q,t) 
ap = — ' 



f d'q D2{q,t) 

J (27r)3 2„q^ ' ^'^^ 



3 15i (27r)3 2iJ,q 

which do not have any finite wavevector singularity. Expressions (23) are obtained, as mentioned before, by truncating 
the perturbation series at the tree level. Extensions to higher orders are straight forward. Illustrative examples of 
higher order diagrams have been shown in Fig. 2. 

Let us now consider various k dependences of Di{k) and D2{k). When the background velocity field is driven by the 

Navier-Stokes equation with a conserved noise (thermal noise) one requires that Di{k) = Dik^, Di = -Di|fc|, giving 
(wj(k, f)Mj(— k, f)) — constant. If we assume similar /c-dependences for (&j(k, 0)&i(— k, 0)) then we require D2{k) ~ 
constant and D2{k) = j^. These choices yield 
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253 f d^q Di\q 

CtD 



aw = — 

3 15 



3 J {2Trf + ^l)q^' 
2834 f d^q e^p[2{aD - aF)\q\t] 



f d- <i exp[2(aj, - aF)\q\t\ , . 

which remain finite even if the system size diverges. 

A fully developed turbulent state, characterised by K41 energy spectra, is generated by Di{k) ^ and Di{k) = 
Dik~'^. In addition if we assume that the initial magnetic fields correlation also have K41 scaling then D2{k) ~ k~^^^ 
and D2{k) = D2k~^/^. If one starts with a K41-type initial correlations for the magnetic fields, then at a later time 
the scale dependence for the magnetic field correlations arc likely to remain same; only the amplitudes grow. Notice 
that the spectra diverge as wavevector A: ^ 0, i.e., as the system size diverges. This is a typical characteristic of fully 
developed turbulence. For such a system we find 



^253 f d^'q Diq-^ _ 2S3 f d^q 
3 7 (27r)3 z/(z/ + 3 J (27r)3 i/' 



2M-i£)i 



.2(1 + P„)g6' 

2^3 4 /• d^q D2ft)g-^/' ^ 253 4 /■ d\ D2(0g-'/' . . 

3 15 J (27r)3 2n 3 15 7 (27r)3 2P„z/ ' ^ ^ 

The notable difference between the expressions Eqs.(24) and (25) for the a coefficients is that the a coefficients diverge 
with the system size if the energy spectra are singular in the infra red limit (as in for fully developed turbulence). 
These divergences are reminiscent of the divergences that appear in critical dynamics [18] which are handled by 
renormalisation group methods. 

In general, at early times (small ap), ctp increases exponentially in time. The growth rate of ap decreases with 
time. Since an and ap have different signs, |(Q!£) + aF)| — > as time t increases. Thus the net growth rate comes down 
to zero. Hence, Eq.(24) and Eq.(25) suggest that the early-time growth and late time saturation of magnetic fields 
take place for different types of background velocity correlations and initial magnetic field correlations. Therefore 
dynamo instability and its saturation are rather intrinsic properties of the 3rfMHD equations with broken reflection 
invariance. One may also note that for i4'41-type of correlations (singular in the infrared limit) one has forward 
cascade of kinetic energy [24]: This is because energy is fed into the system mostly in the large scale (i.e., for small 
k) whereas, dissipation acts primarily in the small scales (large k), resulting into a cascade of energy from the large- 
to small- scales. On the other hand, for correlations smooth in the infra red limit, there is no such cascade. These 
results indicate that the existence of the dynamo mechanism does not require any special background velocity field 
spectrum, though the value of the a-coefficient depends upon it. Our results also suggest that these processes may 
take place for varying magnetic Prandtl number P„i = fi/iy. The above analysis crucially depends on the fact that 
ap and an have opposite signs, which, in turn, imply that Di and D2 have same signs. We have already seen that 
in a physically realisable situation where parity is broken entirely due to the global rotation, £>i and D2 indeed have 
the same sign. 

In the first order smoothing approximation [1,26] in the kinematic limit, to calculate (u x b) one considers only the 
Induction equation as u is supposed to be given. However when one goes beyond the kinematic approximation, one 
has to consider the Navicr-Stokes equation as well. Thus in the first-order smoothing approximation one writes the 
equations for the fluctuations u and b as (to the first order) 

— « V X (u X B) + V X (U X b), (26) 

and 

|^«... + (B.V)b, (27) 

where the ellipsis refer to all other terms in the Navier-Stokes equation and B and U are the large scale {mean field) 
part of the velocity and magnetic fields [1,26]. With these we can write 

d d 

(u X b)i = {cijpUjbp) = {eijpUjBjnj: — Up) + {eijpbpBmj: — bj) = aimBm + ■■■ (28) 

Here the ellipsis refer to non-a terms in the expansion of (u x b) (see Eq.(4)). Thus for isotropic situations a = 
^[— (u.(V X u)) -|- (b.(V X b))] where r is a correlation time. Thus a is proportional to the difference in the fluid 



8 



and magnetic torsalities [4], (fluid hclicity being the same as fluid torsality and magnetic helicity being proportional 
to magnetic torsality) a result obtained in [4,7] using other methods and approximations. Note that Eqs. (21) and 
(23) are very similar to but not exactly the one that were obtained in [7] (in our notations Di is proportional to fluid 
torsality (or fluid helicity) and D2 is proportional to magnetic torsality). Wc ascribe this difference to the essential 
difference between a two-scale approach and our dirgrammatic perturbation theory which, we believe is more suitable 
for handling continuous kinetic and magnetic spectra. 



IV. HYDRODYNAMIC LIMIT OF DYNAMO GROWTH 



We have seen that in Eqs. (25) the a-coefficients diverge in the hydrodynamic (fc — *■ 0) limit which calls for a 
renormalisation group (RG) analysis as a natural extension of our diagrammatic perturbative calculations. In fully 

developed SdMHD, in the steady state, correlation and response functions exhibit dynamical scaling with the dynamic 
exponent 2 = 2/3 [27,28] (for a different approach see [29]), which means renormalised viscosities (kinetic as well as 
magnetic) diverge ~ k~^/^ for a wavenumber k belonging to the inertial range. Even for decaying MHD with initial 
K41-type correlations this turns out to be true [30] where equal time correlations exhibit dynamical scaling with 
z = 2/3. The question is, what it is in the initial transient of dynamo growth [t <C saturation time). We examine this 
in a renormalization group framework. Since we are interested in the early growth, we neglect the Lorentz force and 
work with Eq.(7) inconjunction with the initial magnetic field correlations and noise correlation given by Eq.(ll). As 
before, we assume a statistical steady state for the velocity field. It is well-known that correlations (Mi(k, i)uj(— k, 0)) 
exhibit scaling form k~'''~'^^h{tk^) where x and z are the spatial scaling and dynamical exponents respectively [24] 
where /i is a scaling function. The Galilean invariance of the MHD equations constraints these exponents to obey the 
relation x + z = 1 [24,28,33]. In addition to that, for fully developed turbulence due to non-renormalization of the 
noise-correlators [cf. Eq.(8)] the exponents are fully determined: z = 2/3, x = 1/3, which means the renormalised 
fluid viscosity diverges as fc"^/"^ in the limit wavcvector k ^ 0. During early growth, equal-time magnetic fleld 
correlations (6i(k, t)6j(— k, t)) are expected to exhibit a scaling form k~'^~'^^'"m{tk^) [t <C saturation time) where Xb 
and Zb are the magnetic spatial scaling and dynamical exponents respectively, and to is a scaling function. Similar 
conditions arising from the Galilean invariance and non-renormalization of the initial K41-like magnetic field spectrum 
determines z = Z}j = 2/3 and x = Xft = 1/3- We perform a renormalization group analysis following [19,24,30]. As 
mentioned earlier, the a-term is treated as a perturbation. In a renormalisation-group transformation, one integrates 
ont a shell of modes Ae^' < q < A, and and simultaneously rcscalcs length scales, time intervals and fields through 
X e'x, t e'^i, u — » e'-^u, b e^^h. This has the effect that the nonlinearities are affected only by naive rcscaling 
(this, a consequence of the Galilean invariance of the 3(iMHD equations, essentially implies that the diagrammatic 
corrections to the nonlincarties vanish in the long wavelength limit). The variances Eq.(8), which diverge at low 
wavenumbers remain unrenormalised and thus affected only by rescaling. There are however fluctuations corrections 
to /i and «£) which we evaluate at the lowest order. The resulting RG flow equations for fi and ao, obtained in a 
one-loop calculation are 

t=^[-^-^ + -^ V^(.+V)A^ ]' (2^) 
^ = an[zb - 1 + ^2^^^;^^^], (30) 

where Ai, A2 are numerical constants. Equations (30) and (30) are similar to those presented in Ref. [32] [Eqs. (10.13) 
and (10.14)] but not exactly same. The differences arise mainly (apart from some detail technical differences in the 
perturbation theories involed) from the fact that in Ref. [32] the expressions for the a-coefficients were derived for 
a given variance of the velocity field. In contrast, we use the Navier-Stokes equation, driven by a stochastic force 
of given variance, in place of a given velocity variance. By substituting the value of the exponents in Eqs. (30) and 
(30) we find renormalized (i.e., wavevector dependent) aoik) ~ aDk~^^^, iJ,{k) ~ k~^^^ in the hydrodynamic (fc 0) 
limit. Thus in that limit, the effective dynamo equation takes the form 

^ = {an-,^)k'^% + .... (31) 

where the ellipsis refer to non-linear terms and i refers to the growing mode. Thus, in the hydrodynamic limit, there 
is growth of the magnetic fields only if — > 0. This can happen only if the renormalised magnetic viscosity 
is less than a critical value, set hy an, i.e., the kinetic helicity. In terms of the Ekman number M this condition 
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means M ^ 1/2 for anti dynamo, i.e., no growth, cquivalcntly M "^1/2 for growth of the magnetic fields. This can 
be achieved in two ways, namely by increasing rotation, keeping the magnetic viscosity (or the magnetic Reynolds 
number) constant, or decreasing the magnetic viscosity (i.e., increasing the magentic Reynolds number) for a constant 
rotation. This conclusions are in good agreement with the numerical results of Rcf. [14]. Since renormaliscd magnetic 
viscosity increases with its bare (microscopic) value, it suggests that bare magnetic viscosity must be less than a 
critical value for growth to be possible. Thus our RG results qualitatively explain the numerical results of Kida et 
al [31] who they found that unless magnetic hyperviscosity was less than a critical value there was no growth (it can 
be easily argued that a hypermagnetic viscosity gives rise to a magnetic viscosity in the longer scale and hence their 
result in effect imposes a critical value of the magnetic viscosity). In our model a-effect is proportional to D which 
in turn is proportional to the global rotation frequency. Hence our results suggest that a-effcct is likely to grow with 
increasing rotational speed which is in agreement with the results of Ref. [9] . On the other hand, if the background 
velocity and the initial magnetic field correlators do not have an infra red singularity (i.e., when the correlators ~ k"^) 
there is no fluctuation correction to the magnetic viscosity and to the a-cocfficient resulting in the fact that the 
growth term (oc k) dominates over the dissipation (oc k"^) for sufficiently small wavenumber /e, leading to growth even 
for arbitrarily large magnetic viscosity. Therefore, there is no critical Rm- Thus the effects of the infrared divergences 
that appear in the expressions for the a-cocfficicnts [Eq. (25)] are quite significant: They indicate, as for the driven 
diffusive nonequilibrium systems with diverging kinetic coefficients in the hydrodynamic limit [19,24], divergence of 
time-scales in the hydrodynamic limit. Since, the a-term in Eq. (7) is proportional to wavenumber fc, the time-scale 
of growth of the mode with wavenumber k is 0{ak). This remains true, even in the hydrodynamic limit, for the 
case when there is no divergence in the a-coefhcients. In contrast, when the a-coefffcient diverge in the infra red 
limit, the growth rate changes qualitatively from its linear dependence on wavenumber k in the hydrodynamic limit. 
For example, with the the background velocity correlations and the initial magnetic field correlations given by Eq. 
(8), the a coefficients diverge as A;"^/"^ in the long wavelength limit. Hence, the effective growth rate is changed to 
a{k)k ~ fc^/^. A full self-consistent calculation (when feedback due to the Lorentz force cannot be neglected) for the 
a-coefficients require simultaneous solutions of the self-consistent expressions for magnetic Prandtl number, magnetic- 
to kinetic- energy ratio and the a-coefficients which can be handled in our scheme of calculations. The self-consistent 
solutions are expected to be influenced by the degree of crosscorrelations between the velocity and magnetic fields 
[33]. 

So far, we have assumed that both Di{k) and Di{k) have the same infra red singularity {Di{k) ~ k and 
Di{k) ~ k~^/^). This need not be the case always. However, if Di{k) is non-singular then does not diverge. 
As a result, the growth rate is just ajjk even in the hydrodynamic (long wavelength) limit. Effective dissipation, 
however, will still be ~ k'^^^ and thus it will dominate over 0{k) growth. Therefore, there will be no growth in the 
hydrodynamic limit. Thus, our analyses suggest that in any fully developed turbulent system with a-effect, helicity 
spectrum (given by Di{k)) should be as singular as the kinetic energy spectrum (given by Di{k)). 

V. CONCLUSIONS 

In conclusions, we have calculated expressions for the a-coefficients in a diagrammatic perturbation theory on 
a minimal model for arbitrary background velocity and initial magnetic field correlations, and fluid and magnetic 
viscosities. We show that the parity breaking parts of the velocity and magnetic field variances must have the same 
sign, which is the case in any physical system. Wo explicitly show that the processes of early growth and late-time 
saturations may take place independent of any special velocity and initial magnetic field correlations. Even though 
our explicit calculations were done by using simple initial conditions for the calculational convenience, the results 
that wc obtain are general enough and it is apparent that the feedback mechanism is qualitatively independent of the 
details of the initial conditions and force correlations, one may note that for one of the force/initial correlations there 
is no kinetic energy cascade in the conventional sense but we still find dynamo action. It is quite reasonable to expect 
that our results should be valid for more realistic initial conditions also. In effect we have explicitly demonstrated 
the robustness and generality of the dynamo mechanism and that the dynamo mechanism is an intrinsic property 
of the 3dMHD equations. We have also shown, within our RG analysis, that the magnetic viscosity should be less 
than a critical value for growth of magnetic fields a result which was previouly observed in numerical simulations. Wc 
conclude the existence of a critical Ekman number for K41-type correlations: We find growth only when M < 1/2, 
confirming recent numerical results. This is easily understood in our framework. The issue of divergent effective 
viscosities in the inertial range assumes importance as it may help to overcome some of the non-linear restrictions as 
discussed by Vainshtein and Cattaneo [6]. A system of magnetohydrodynamic turbulence in a rotating frame, after 
the saturation time (i.e., after which there is no net growth of the magnetic fields) belongs to the universality class 
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of usual three-dimensional magnetohydrodynamic turbulence in a laboratory. This can be seen easily in both the lab 
and the rotating frames; the critical exponents characterising the correlation functions can be calculated exactly by 
using the Galilean invariance and noise- nonrenormalisation conditions [24,28]. An important question, which remains 
open for further investigations, is the multiscaling properties of the velocity and the magnetic field structure functions 
at various stages of the growth of the magnetic fields. In what concerns an experimental observation of our results, 
one should add that even though it is not easy to verify our results in an experimental set up, numerical simulations 
of Eqs.(6) and (7) with the variances (11) with different fc-dependences can be performed to check these results. 
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